如何用stata计算基尼系数?

您所在的位置:网站首页 stata beta命令 如何用stata计算基尼系数?

如何用stata计算基尼系数?

#如何用stata计算基尼系数?| 来源: 网络整理| 查看: 265

Source: 本文根据如下论文翻译整理而得,特此致谢! Owen O'Donnell, Stephen O'Neill, Tom Van Ourti, Brendan Walsh, 2016, Conindex: Estimation of Concentration Indices, Stata Journal, 16(1): 112–138. 1. 背景介绍

「集中度指数」 衡量的是一个变量相对于另一个变量的分布的不平等 (Kakwani, 1977) ,为评估不平等,标准集中度指数与广义集中度指数等多种集中度指标已被提出。

本推文将要介绍的conindex命令提供了一种简单的统一方法来估算各种集中度指数及其标准误差。 它可以测量基本变量的横截面不平等,也可以比较重复的横截面或面板数据各个时期的不平等,还可以用它来估计单变量不平等的等级相关指数,例如「基尼系数」和「广义基尼系数」。

与其他计算等级相关的不平等指数的命令,如 concindc (Chen, 2007) 和 glcurve (van Kerm 和 Jenkins, 2007) 等相比, conindex 命令估计了一系列集中度指标,这使分析人员可以选择一个适合所关注变量的测量属性的指标,并且该指标应符合其关于不平等性的规范性原则。

2. 集中度指数2.1 理论部分

conindex 命令可以计算的各种不平等指数, 为了简化说明,将兴趣变量称为健康,将排名变量称为收入,则与收入相关的健康不平等状况可以通过绘制从最贫穷到最富有的个人的健康累积比例来进行评估。具体的集中度指数测度定义如下:

C(h \mid y)=\frac{2 \operatorname{cov}\left(h_{i}, R_{i}\right)}{\bar{h}}=\frac{1}{n} \sum_{i=1}^{n}\left\{\frac{h_{i}}{\bar{h}}\left(2 R_{i}-1\right)\right\} \\

其中,

h_{i} 是用于衡量不平等的变量,例如,健康。C 的范围从 (1-n) / n ,即最大的“亲贫” (pro-poor) 不平等 (也就是说, 所有健康都集中在最贫穷的个体上) 到 (n-1) / n ,即最大的“亲富” (pro-rich) 不平等。 公式 (1) 显示,集中度指数可以解释为 (建康) 份额的加权平均值,其权重取决于分数 (收入) 等级 \left(2 R_{i}-1\right) 。 当在用于排名的同一变量中测量不平等时,单变量不平等的基尼系数度量作为集中指数的一种特殊情况出现。以上标准集 中指数衡量相对不平等,并且对于关注变量 (建康) 的等比例变化是不变的。但在另一个极端,绝对不变性对应于不平等度 量,当两个分布显示相对不平等的水平相同时,均值越高的绝对不平等就越大。广义集中度指数 GC 可以表示为:

\operatorname{GC}(h \mid y)=\frac{1}{n} \sum_{i=1}^{n}\left\{h_{i}\left(2 R_{i}-1\right)\right\} \\

范围介于 \bar{h}\{(1-n) / n\} (最亲贫) 和 \bar{h}\{(n-1) / n\} (最亲富) 之间。

在测量水平允许对感兴趣变量进行转换 (标量,序数,基数,比率或固定比例) 的情况下,标准集中度指数和广义集中度指数不一定是不变的或相等的,还可以演化出:

广义 (基尼) 集中度指数Wagstaff 和 Erreygers 归一化集中度指数……就在这里就不一一详细介绍,感兴趣的可以参考 Owen O’Donnell et al. (2016) 的文章。

2.2 推断与估计 实际上,每个与等级相关的不平等指数都可以表示为关注变量 \left(h_{i}\right) 与排序变量的分数等级 \left(R_{i}\right) 之间的协方差的转换。例 如,标准集中度指数等于相关变量的均值除以协方差的两倍[见公式(1)]。而简单最小二乘回归的斜率系数是协方差除以回 归元的方差,所以每个不平等指标都可以通过对目标变量在等级上的变换进行回归得到。例如,标准集中度指数是对如下 回归模型中系数 \alpha_{1} 的最小二乘估计:

\frac{2 \sigma_{R}^{2}}{\bar{h}} h_{i}=\alpha_{0}+\alpha_{1} R_{i}+\epsilon_{i} \\

其中, \sigma_{R}^{2} 是 R 的方差,而 \epsilon_{i} 是误差项。 \alpha_{1} 的最小二乘估计标准误差用作集中度指数的估计标准误差。

3. conindex 命令语法与实例3.1 命令语法

conindex 的第一版本发布于 2016 年,目前更新版为 2018 版本。输入 findit conindex 可以搜索到该命令的最新版本,找到最新版本后,点击 (click here to install) 即可。亦可使用命令来安装:

. net install st0427_2.pkg

完成安装后,输入 help conindex 即可查看其完整帮助文件。

conindex 的语法语法结构如下:

conindex [varlist] [if] [in] [fweight aweight pweight], RANKvar(varname) [CLUSter(varname)] [robust] [truezero] [generalized] [v(#1)] [beta(#1)] [bounded] [LIMits(#1 #2)] [WAGstaff] [ERReygers] [graph] [loud] [COMPare(varname)] [KEEPrank(string)] [svy]

没有默认要计算的指数,要通过选项定义要计算的指数,例如,广义集中度指数是通过同时指定 truezero 选项和 generalized 选项而得出的,其他指数的具体指定说明可参考 Owen O’Donnell 等 (2016) 的文章,也可以通过在 stata 命令窗口输入 help conindex 命令得到。 同时,可以使用 graph 选项获得 (广义的) Lorenz 和 (广义的) 集中度曲线,也可以使用 xtitle() 和 ytitle() 选项设定横坐标和纵坐标的标题。

3.2 Stata 实例

由于 Owen O’Donnell 等 (2016) 所使用的示例数据无法获取,下面我们使用模拟数据说明 conindex 的功能,具体命令如下:

Note: 下面的模拟过程中,未设定种子值。因此,你自己测试的时候,结果会与我们呈现的结果略有差异。 . clear . set obs 1000 //设定1000个观测值 //生成两个模拟变量 . gen rank = _n . gen Y1 = invnormal(uniform())^2 //标准集中度估计 . conindex Y1, truezero rankvar(rank) -------------------------------------------------------------+ Index: | No. of obs. | Index value | Std. error | p-value | --------+-------------+-------------+-------------+----------| CI | 1000 | .04350093 |.02544088 | 0.0876 | -------------------------------------------------------------+ //广义集中度估计 . conindex Y1, generalized truezero rankvar(rank) -------------------------------------------------------------+ Index: | No. of obs. | Index value | Std. error | p-value | --------+-------------+-------------+-------------+----------| Gen. CI | 1000 | .04519981 |.02643444 | 0.0876 | -------------------------------------------------------------+ //生成受限变量 ,也就是有上下限的变量 . gen BinaryY=invnormal(uniform())>0 . sum BinaryY Variable | Obs Mean Std. Dev. Min Max ----------+-------------------------------------- BinaryY | 1,000 .466 .4990923 0 1 //此时,需要指定bounded和limits选项 . conindex BinaryY, bounded limits(0 1) rankvar(rank) --------------------------------------------------------------+ Index: | No. of obs. | Index value | Std. error | p-value | ----------+-------------+-------------+------------+----------| CI | 1000 | -.01971674 |.01955377 | 0.3135 | --------------------------------------------------------------+ //也可以通过指定wagstaff和erreygers选项来获得归一化浓度指数 . conindex BinaryY, bounded limits(0 1) wagstaff rankvar(rank) -----------------------------------------------------------------------+ Index: | No. of obs. | Index value | Std. error | p-value | -------------------+-------------+-------------+------------+----------| Wagstaff norm. CI | 1000 | -.03692273 |.03661754 | 0.3135 | -----------------------------------------------------------------------+ . conindex BinaryY, bounded limits(0 1) erreygers rankvar(rank) -----------------------------------------------------------------------+ Index: | No. of obs. | Index value | Std. error | p-value | -------------------+-------------+-------------+------------+----------| Erreygers norm. CI | 1000 | -.036752 |.03644822 | 0.3135 | -----------------------------------------------------------------------+

详细内容参见连享会推文

专题:Stata教程 Stata:各类集中度指数估算-广义基尼Gini系数相关推文 Note:产生如下推文列表的 Stata 命令为:. lianxh 系数 . songbl 系数安装最新版 lianxh/ songbl 命令:. ssc install lianxh, replace 专题:Stata教程 Stata:自己动手做组间系数差异检验-bootstrap-bdiffStata:各类集中度指数估算-广义基尼Gini系数专题:Stata命令 敏感性分析B-Stata实操:控制变量内生时的系数敏感性分析-regsensitivity敏感性分析A-理论基础:控制变量内生时的系数敏感性分析-regsensitivityStata:系数稳定性分析-psacalc专题:Stata绘图 Stata绘图:回归系数可视化-multicoefplotforest-森林图:分组回归系数可视化Stata绘图:回归系数可视化-论文更出彩Stata绘图:随机推断中的系数可视化专题:结果输出 Stata结果输出:pwcorr_a输出相关系数矩阵pwcorr_a:输出相关系数矩阵至 Word 和 ExcelStata:回归结果中不报告行业虚拟变量的系数专题:回归分析 Stata:分组回归系数比较的新思路0.0005:估计系数太小怎么办?Stata:系数为何不显著?GIF-演示-OLS-的性质.md正确姿势:回归系数的解释与评估Stata: 实时估计个股贝塔(beta)系数多元回归系数:我们都解释错了?图示线性回归系数:Frisch-Waugh定理与部分回归图Stata: 如何检验分组回归后的组间系数差异?如何比较解释变量的系数相对大小?Stata: 获取分组回归系数的三种方式专题:面板数据 异质性分析:系数平滑可变模型面板变系数模型:每家公司都有一个斜率专题:IV-GMM IV-工具变量法:第一阶段系数符号确定时的小样本无偏估计专题:交乘项-调节 我为什么越跑越肥?交乘项系数的直观解释专题:时间序列 DCC-MGARCH:动态条件相关系数模型 (R+Stata)CCC-MGARCH:常值条件相关系数模型 (R+Stata)专题:风险管理 Stata:普通年金现值系数表


【本文地址】


今日新闻


推荐新闻


CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3